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Supersonic  Nozzle  Design.  (Project  NOL  159) 


Abstract:  An  analytic  methcd  for  designing  nozzles  is  described  herein. 

This  method  (Nil son' 3)  allows  for  the  nonuniformity  of  flow 
across  the  throat. 


Foreword:  This  supplement  covers  material  that  was  the  subject  of  a 

course  on  nozzle  design  given  by  the  author  and  Dr.  7f.  H.  Heybey. 
The  data  and  conclusions  presented  here  are  for  the  use  of  the 
personnel  of  the  Naval  Ordnance  Laboratory.  They  do  not  neces- 
sarily represent  the  final  opinion  of  the  Laboratory. 
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Baltimore,  273  (1930)7~ 

Enel: 

(A) 

Figure  1 
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Table  I 
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I.  Nilson1 s Method 
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<*' 


1 
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1.  Various  people  have  developed  exact  analytic  methods  for  design- 
ing nozzles;  however,  the  author  is  only  sufficiently  familiar  with  that 
of  Nilson  (references  a and  b)  to  discuss  it  herein. 

2.  Essentially  his  method  consists  of  finding  an  exact  solution  for 
flow  in  the  region  of  the  throat  and  then  converting  the  divergent  super- 
sonic flow  thus  produced  into  a uniform  stream  at  the  desired  Mach  number. 
This  conversion  into  a uniform  stream  is  similar  to  the  ccrr/'ision  of 
source  flow  into  a uniform  stream  in  the  Foelsch  mu  ..hod;  the  difference 
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between  the  two  methods  being  that  Foelsch  assumes  that  he  will  have 
source  flow  while  Nilson  has  an  exact  solution  and  knows  what  velocity 
distribution  his  contour  will  produce.  Both  consider  only  symmetrical 
nozzles. 


3.  More  specifically,  Nilson  assumes  a given  velocity  distribution 
along  the  nozzle  centerline.  He  expresses  his  ordinates  and  velocity  as 
infinite  series,  the  coefficients  of  which  he  determines  from  the  dynamical 
equations.  This  method  of  designing  nozzles  is  naturally  limited  to  cases 
in  which  the  series  converge  within  the  desired  accuracy;  (i.e.,  at  more 
or  less  low  Mach  numbers).  To  insure  the  proper  convergence  the  present 
author  evaluated  the  coefficients  of  the  third  terms;  Nilson  gave  only 
the  first  two. 


4.  Let  us  adopt  the  following  notation: 
x,  y <=  rectangular  coordinates,  X axis  coincident  with  nozzle  centerline 
r.Q  - magnitude  and  direction  (relative  to  X axis)  of  velocity  vector. 


^ H » coordinate  system  employing  equipotential  line  and  streamlines 
' (see  reference  a). 


M 


-j / 
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where  o*  is  the  critical  sonic  speed 


ln=  £ = 

A*  M 


Hi 

2. 


A bar  over  a variable  indicates  that  the  particular  variable  is  taken  on 
the  nozzle  centerline. 


5.  The  solution  is  assumed  in  terms  of  the  following  infinite 
series : 


* - F v-  yj  V x¥  y * 


Q-  Q,^-t  d3 


a) 

(2) 

(3) 

(4) 
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The  coefficients  were  found  to  be 


(prime  denotes  differentia- 
tion with  respect  to  |*  ) 


(5) 

(6) 

(7) 

(8) 


(9) 


6.  The  function  In  ~ occurring  in  the  preceding  formulae 

is  the  assumed  centerline  velocity  distribution.  For  a distribution 
A - » an  approximately  hyperbolic  contour  is  generated.  If 

this  flow  field  is  now  cut  off  along  any  Mach  line  and  the  flow  con- 
cried  into  a uniform  and  parallel  stream,  t are  will  be  a discontinuity 
in  streamline  curvature  at  the  junction  of  tlx.se  two  flows.  Such  a 
discontinuity  is  not  desirable;  but,  however,  what  distribution  must 
be  chosen  to  yield  streamlines  possessing  inflexion  points  is  not  known. 
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For  low  Mach  numbers,  at  least,  the  change  in  curvature  is  not  too 
serious  providing  proper  discretion  is  used.  How  much  curvature  discon- 
tinuity can  be  safely  used  remains  to  be  determined  from  experiment. 

7.  Using  h~  , we  can  simplify  some  of  the  coefficients 

in  the  preceding  equations.  By  definition 

h=  -1  (for  y = 1.4) 

2/6  fv\ 

and  solving  this  equation  from  K in  terms  of  h,  we  obtain 

which  for  h ~ /-/  J"  becomes 
then,  approximately 

(OJH7I&4  03  /Z.JIZQ  f-L'")  (6«) 


(uf'p.ZWOOi  O.M£b$i  bVi7lcj>-) 
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8.  A sketch  of  the  nozzle  illustrating  notation  is  given  on 
Figure  1.  The  design  of  a nozzle  is  started  at  pt  L*  , the  point  on 
the  nozzle  axis  at  which  the  final  (test  section)  Mach  number  is  first 
achieved.  Let  this  Mach  number  be  . The  corresponding  area 

ratio  can  then  be  easily  obtained.  The  coordinates  of  this 

point  are  seen  to  be  X-^r  ^ where,  for  h~/-f-ez  , //^  — / 1 

and  . Using  these  values  as  initial  conditions,  we  can 

compute  the  Mach  line  which  is  given  by  the  equation 


il 


. m 


(14) 
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9.  This  equation  must  be  integrated  numerically.  The  Kutta-Runge 
method,  described  in  reference  (c),  can  be  used.  Equation  14  is  in  the 
form 
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One  integrates  step  by  step.  First  compute 

k=f(i - f?,  y |v^| 

4=^/^  4/, 

4=^y  ^ 7, 


(14a) 


(15a) 

(15b) 

(15c) 

(15d) 


Then 


=£{%■* 2 4+Z-4+4J  ' 


(15e) 


as  a function  of  /^)"/  can  be  obtained  from  Table  I and 
41*  from  equation  3. 

/ft* 

10.  The  integration  is  extended  to  a suitable  point  as  previously- 
discussed.  Holding  the  value  of  h at  thiB  point  constant,  we  can 
obtain  from  equations  1 and  2 the  contour  of  the  nozzle,  a streamline, 
upstream  of  this  point.  The  contour  shculd  be  extended  approximately 
again  as  far  as  the  minimum  section.  The  coordinates  of  the  minimum 
section  can  be  found  from  the  relation 

%“) 

- - ' ' q. 

One  should  note  that  for  1-0  / x = -am'jib  rj  . 


(16) 


11.  The  coordinates  V, y of  the  remaining  (terminal)  portion 
of  the  nozzle  can  be  determined  from  the  following  equations  (see 
Figure  1) : 


+ l cos  (0, 
y ~ y,  i l Sin  ( 9/ 
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M as  a Function  of\|b 


M 

M 

4E 

\ * 

& 

JL 

.00 

1.00000 

.43 

1.51276 

.87 

2.04793 

.00 

1.0000c 

.01 

1.01099 

.44 

1.52519 

.88 

2.0:5959 

-.01 

.9890 L 

.02 

1.02204 

.45 

1.53761 

.89 

2.07122 

-.02 

.97822 

.03 

1.03316 

.46 

1.55004 

.90 

2.08283 

-.03 

.96744 

.04 

1.04434 

.47 

1.56246 

.91 

2.09440 

-.04 

.95673 

.05 

1.05558 

.48 

3.57488 

.92 

2.10595 

-.05 

.94609 

.06 

1.06688 

.49 

1.58730 

.93 

2.11746 

-.06 

.93552 

.07 

1.07824 

.50 

1.59971 

.94 

2. ‘42894 

-.07 

.92502 

.08 

1.08965 

.51 

1.61212 

.95 

2.14040 

-.08 

.91461 

.09 

1.10112 

.52 

1.62452 

.96 

2.15182 

-.09 

.90426 

.10 

1.11265 

.53 

1.63691 

.97 

2.46321 

-.10 

.89398 

.11 

1.12422 

.54 

1.64929 

.98 

2.17457 

-.11 

.88381 

.12 

1.13585 

.55 

1.66166 

.99 

2.18590 

-.12 

.87370 

.13 

1.14752 

.56 

1.67403 

1.00 

2.49720 

-.13 

.86367 

.14 

1.15925 

.57 

1.68638 

1.01 

2.20347 

-.14 

.85372 

.15 

1.17102 

.58 

1.69872 

1.02 

2.21970 

-.15 

.84386 

.16 

1.18283 

.59 

1.71104 

1.03 

2.23091 

-.16 

.83408 

.17 

1.19469 

.60 

1.72336 

1.04 

2.24208 

-.17 

82433 

.18 

1.20658 

.61 

1.73565 

1.05 

2.25322 

-.18 

.81476 

.19 

1.21852 

.62 

1.74793 

1.06 

2.36433 

-.19 

805  23 

.20 

1.23050 

.63 

1.76020 

1.07 

2. 27541 

-.20 

7957? 

.21 

1.24251 

.64 

1.77244 

1.08 

2 23645 

-.21 

.7864 2 

.22 

1.25455 

.65 

1.78467 

1.09 

2.  2)747 

-.22 

.77715 

.23 

1.26663 

.66 

1.79688 

1.10 

2.  345 

-.23 

.76796 

.24 

1.27874 

.67 

1.80907 

1.11 

2.  U940 

-.24 

.75886 

.25 

1.29008 

.68 

1.82124 

1.12 

2.  <3032 

-.25 

.74985 

.26 

1.30305 

.69 

1.83339 

1.13 

2. 34120 

-.26 

.74093 

.27 

1.31524 

.70 

1.84552 

1.14 

2.  3 06 

-.27 

.73209 

.28 

1.32746 

.71 

1.85763 

1.15 

2. 26288 

-.28 

.72334 

.29 

1.33969 

.72 

1.86971 

1.16 

2. 

-.29 

.71469 

.30 

1.35197 

.73 

1.88177 

1.17 

2. 3?443 

-.30 

.70612 

.31 

1.36426 

.74 

1.89380 

1.18 

2.39516 

.32 

1.37656 

.75 

1.90502 

1.19 

2. 40585 

.33 

1.38889 

.76 

1.91780 

1.20 

2.41651 

.34 

1.40123 

.77 

1.92977 

1.21 

2.42715 

.35 

1.41358 

.78 

1.94170 

1.22 

2.43725 

.36 

1.42595 

.79 

1.95362 

1.23 

2<  44831 

.37 

1.43832 

.80 

1.96550 

1.24 

2.45885 

.38 

1.45071 

.81 

1.97736 

1.25 

2. 46936 

.39 

1.46311 

.82 

1.98919 

1.26 

2. 47983 

.40 

1.47552 

.83 

2.00099 

1.27 

2. 49027 

.a 

1.48793 

.84 

2.01277 

1.28 

2,30068 

.42 

1.50034 

.85 

2.02452 

1.29 

2.  51106 

.43 

1. 51276 

.86 

2.03624 

1.30 

2.52ia 

.87 

2.04793 

*Thie  Table  appeared  in  reference  (a). 
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